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Abstract. We show that a model for radio source dynam- 
ics we had earlier proposed can readily reproduce the rela- 
tionship between the radio power division separating the 
two Fanaroff-Riley classes of extragalactic radio sources 
and the optical luminosity of the host galaxy, as found 
by Owen and Ledlow (1994). In our scenario, when less 
powerful jets eventually slow down to the point that the 
advance of the working surface (i.e., hotspot) becomes 
subsonic with respect to the external gas, the jet's col- 
limation is severely weakened. This criterion distinguishes 
the powerful and well collimated FR II sources from the 
weaker sources producing the less collimated FR I type 
morphologies. 
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1. Introduction 

A landmark in the study of extragalactic radi o sour ces 
was the demonstration by Fanaroff and Riley (1974) of 
the existence of a relatively sharp morphological transition 
at a radio luminosity of corresponding to Pi78MHz — 
2.5 x lO 26 /^ 2 W Hz -1 . The great majority of sources be- 
low this luminosity (FR I type) are characterized by hav- 
ing diffuse radio lobes, with their brightest regions within 
the inner half of the radio source, and so can be called 
edge-dimmed. On the contrary, more powerful sources are 
usually straighter, exhibit edge-brightened (FR II) mor- 
phology, and typically contain hotspots near the outer 
edges of their radio lobes. 

More recently it was realized that the critical radio lu- 
minosity separating the FRI and FR II actually increases 
with the optical luminosity of the host elliptical galaxy, 
so that Pft oc io P 6 t 5 (as measured from Fig. 1 of Ledlow & 



Owe n |1 996| ); also see 0wen & White ( |l99l[ ); Owen & Led- 
low fll994j ) for earlier indications of this effect. Although 
many detailed differences between the host galaxies of the 
FR I and FR II sources have been discovered (e.g., Baum 
et al. 1995; Zirbel 1997; for a recent summary, see Gopal- 
Krishna & Wiita [2000 ), it is fair to state that a key fact 
is that the more luminous the host galaxy is, the more 
powerful the radio source must be in order to attain the 
FR II morphological classification. 

The most widely known explanation for the difference 
between the FR I and FR II sources is that, while the 
jets in both cases start out moving at very high (rela- 
tivistic) speeds, those in FR II sources remain that way 
out to multi-kpc distances, while those in the FR Fs de- 
celerate to much lower speeds within a few kpc of the 



Komissarov 



1994 



Laing 



core (e.g., B egelm an 1982; Bickncll 1984 ; De Young 



1992 



1996 



Laing et al. |1999J) 



1993 



Detailed models for deceleratin g relativistic jets were de- 
veloped by Bicknell ( 1994 , 1995 ). The core of his argument 
is that relativistic jets, once they come into pressure equi- 
librium with the external interstellar medium (ISM) or 
intracluster medium (ICM), will become strongly Kelvin- 
Helmholtz unstable and then entrain substantial amounts 



of cold ambient material. Bicknell ( 1995 ) showed that this 
instability typically set in as the internal jet Mach number 
approaches ~ 2, corresponding to internal bulk velocities 
having dropped to about 0.6c. 

Of course, many other explanations for the FR I/FR 



II dichotomy have been proposed (e.g., Baum et al. 1995 
Reynolds et al. |l996a| , |1996b| ; Meier |1999| : Valtoncn 



Heinamaki 2000; see the recent discussion by Gopal- 



Krishna & Wiita 2000). Some of these proposals involve 
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fundamental differences in the nature of the jets (elec- 
tron/positron vs. electron/proton plasma) or of the cen- 
tral engine (black hole spin, type of accretion disk). The 
viability of any member of this class of explanations is 
challenged by the existence of HYbrid MOrphology Radio 
Sources, or HYMORS, which have a clear FR I morphol- 
ogy on one side of the host galaxy, but a distinct FR II 
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morphology on the other side (Gopal-Krishna & Wiita 

( poooj , pooi| ). 

An alternative approach posits that FR I and FR II 
sources differ primarily in the importance of the beam 
thrust relative to the basic paramete rs of the ambi- 
ent medium ( Gopal -Krishna & Wiita |l988| (GK W88); 



Gopal-Krishna |1991| (GK91); Gopal-Krishna et al. |1996 



(GKWH96); Blandford 1996). In this version of the decel- 
eration scenario, the emphasis is on the slowing of the ad- 
vance of the hotspot, or working surface, at the end of the 
jet, rather than on the slowing of the bulk flow within the 
jet (e.g., Bicknell 199J). When the hotspot's advance be- 
comes transonic relative to the ambient medium, its Mach 
disk weakens considerably due to the fall in ram pressure, 
and the jet becomes decollimated; this soon leads to an 
FR I morphology. Further expansion of such jets beyond 
this point is expected to be in the form of a plume, as dis- 
cussed in GKW88, where it was also shown that for rea- 
sonable values of jet and galaxy parameters (density, core 
radius and temperature of the ISM) such a jet flaring can 
occur within several kpc of the core for jet powers below 
about 10 43 erg s _1 , consistent with the observations of the 
FR I/FR II division (GKW88; GK91; GKWH96). Also, in 
this picture, the concomitant dimunition/cessation of the 
'backflow' of the beam plasma would cause a depletion 
of the protective sheath of beam plasma around the jet, 
facilitating entrainment of the ISM material into the jet 
flow. 

Bicknell ( |1995 ) expanded his model to account for the 
Owen-Ledlow transition in the PR-L opt plane. A key in- 
gredient in his scenario involved connecting the jet dynam- 
ics to empirically established relations between the optical 
magnitude of the galaxy on the one hand, and the soft 
X-ray luminosity, core radius, and central velocity disper- 
sion of elliptical galaxies on the other hand. In the present 
study, we follow Bicknell and adopt these same empirical 
relationships. Bicknell then ties the galactic parameters 
to those of the jet by demanding that the establishment 
of a pressure balance between the jet and the external 
medium corresponds closely to the location of the jet's in- 
ternal transonic transition. By doing so, he was able to 
derive a formula connecting beam power to the optical 
luminosity of the host galaxy. In this paper we follow a 
different approach in this latter stage. 

Additional assumptions on the efficiency of the con- 
version of the jet's energy into total radio luminosity 
(> 1 — 2%), and thence into monochromatic radio emis- 
sion, allowed Bicknell ( 1995| ) to obtain a fairly good fit 
to the slope of the Owen-Ledlow division, finding that 
P£ oc iopt) as wen as getting a decent fit to the inter- 
cept; however, there were quite a few poorly constrained 



parameters in his model. Nevertheless, Bicknell (1995) ar 



gued that the the slope of the relation was rather insensi- 
tive to the likely uncertainties in parameters. He further 
argued that the intercept would tend to be driven towards 
better agreement with the data for plausible variations in 



those parameters, which included: the ratio of the jet's life- 
time to the synchrotron cooling time of the highest energy 
electrons; the ratio of the upper and lower cutoff energies 
for the electron distribution; the ratio of jet radius to its 
length at the transition distance; the ratio of the product 
of the jet pressure and square of its radius evaluated at 
the transition radius to that of the external medium eval- 
uated at its core radius; the bulk velocity of the jet at the 
transition point; the low frequency radio spectral index; 
and the value of an integral (discussed below) which de- 
pends on an upper cut-off radius for the X-ray emitting 
halo of the host galaxy. 



The new magnetic switch model (Meier 1999 ) funda- 
mentally distinguishes FR I from FR II sources through 
the different speeds of rotation of the magnetic field lines 
of their central engines, which are tied to the different 
spin rates of their supermassive black holes. While this 
scenario can also produce a slope for the — L opt di 



viding line close to that found by Ledlow & Owen (1996) 
(Meier gggf ) we recall that it is hard to reconcile this 
scheme with the existence of HYMORS (Gopal-Krishna 
k Wiita |200C| , |200l| ). Vent uring beyond the standard nu- 
clear jet paradigm, gravitational slingshot models can in 
principle be compatible with the Owen-Ledlow relation 
(Valtonen & Heinamaki 2000| ), and could also account for 
HYMORS (Gopal-Krishna & Wiita pOOg), 



2. The Model 

In an earlier study involving "weak headed quasars", 
where prominent one-sided jets are not seen to terminate 
in a conspicuous hotspot, we have argued that dissipa- 
tion of jet power, as s uggested by Swarup et al. ( 1982 ) 
and Saikia et al. (1983), may actually not be responsible 
for the lack of terminal hotspots (GKWH96). Taking a 
clue from the fact that no two-sided jets are seen in these 
sources, we argued that the lack of a hotspot could be best 
explained through the onset of a jet's decollimation when 
the hotspot's (or, nearly equivalently, the bow shock's) ve- 
locity becomes transonic relative to the external medium 
(GKWH96). This successful explanation of this class of 
radio sources motivates us to seek an explanation for the 
Owen-Ledlow relation in terms of a similar scenario. This 
model makes no assumptions about the relativistic nature 
of the bulk velocity of the internal jet fluid, which may be 
gradually decelerating as the head of the jet advances. 

Our study makes use of the same empirical relations 
between the elliptical's blue magnitude, Mb, and i ts sof t 
X-ray emission, Lx (Donnelly, Faber and O'C onncll 1990 ), 
stellar velocity dispersion, a (Terlevich et al. [1981 — the 
Faber-Jackson relation), and X-ray core radius, a (Kor- 
mendy 1987 ), as employed by Bicknell ( |1995|) , assuming 
H = 75 km s" 1 A/r ~~ " : 



Mpc" 



log L x = 22.3 - 0.872Mb, 
log a = 5.412 - 0.0959Mb, 



(1) 
(2) 
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log a = 11.7-0.436A/b- 



(3) 



However, we adopt the initially conical jet mod el based 
on our earlier work (Gopal-Krishna & W iita [1987 ), which 
was generalized in Gopal-Krishna et al. ( 1989 , GKWS89) 
to allow for bulk relativistic jet flow whereby the hotspot 
motion can also be relativistic. Ram pressure balance, 
under these conditions, gives the hotspot velocity, v, as 
a function of distance, D, from the central engine, as 
(GKWS89): 



v(D) 



Xc[l + (D/a 



.215/2 



D + X[l + (D/a) 2 } 5 / 2 ' 



(4) 



Here X — (4Li,/TTC 3 9 2 n m P n) 1 ^ 2 , where is the jet 
(beam) power, 9 is the jet's effective opening angle, m p 
and \x have their usual definitions, and uq is the central 
density of the ISM, which is taken to fall off as (e.g., For- 



man et al. 1985; Canizares et al. 1987 Conway 2001) 



n(D) = 



n 



[1 + (D/a) 



215 • 



(5) 



For our calculations we choose 9 = 0.1 rad, which is cer- 
tainly justified as appropriate for the inner jet regions, 
particularly since it is now recognized that the jet's thrust 
acts on a larger area than the instantaneous hot spot (e.g., 
Scheuer's 1982] 'dentist-drill', as supported by many three- 
dimensional numerical simulations: Norman 1996: Clarke 



1996| ; Kaiser & A lexan der |1997] Hooda & Wiita |1998j ). Fol- 
lowing Bicknell ( 1995| ), we further assume that the X-ray 



emitting ISM gas temperature is tied to the central stel- 
lar velocity dispersion, er, via, kT = 2fl[im p a 2 / '5, which 
is Fall's ( |1987[ ) relation between circular velocity and a. 
Here S is roughly 0.75, as inferred from s oft X- ray images 
of ne arby elliptical s (e.g ., Forman et al. 1985 ; Canizares 
et al. 1987| ; Sarazin 1986 ), as well as dyna mical m odels of 
radio sources (e.g., Gopal-Krishna & Wiita 1991 ; Conway 
200ll ). 

Studies of many sub-classes of young radio sources, in- 
cluding Compact Symmetric Objects (e.g., Owsianik & 



Conway [1998 ; Conway 2001 ), Gigahertz Peaked Sources 



(Carvalho 1999) and 



sources (Gopal-Krishna & Wiita 1991 



Compact Steep Spectrum radio 
Jeyakumar et al. 



2001) have indicated that the mean density of the gas in- 
teracting with the radio lobes in the inner ~ 1 kpc or so 
is typically a few atoms cm -3 , or several times that esti- 
mated from X-ray emission. This evidence in favor of a sig- 
nificant contribution of cooler gas to the ambient medium 
(at least in the inner portions of elliptical galaxies) is fur- 
ther supported by studies of the linear-size distribution 
of these radio sources (O'Dea & Baum 1997). We shall 
parameterize the relation between the confining central 
density, no, used in Eq. (5) and that derived from X-ray 
measurements, nx as follows: no — nnx, with k > 1; the 
fiducial value we use below is k = 3. 

The radio data cited above as well as the relation be- 
tween Mb and the X-ray core radius (Eq. 3) lead to values 



of of a ~ 1 kpc. Therefore the ISM density begins to ap- 
proach the large scale behavior, and thus declines quite 
rapidly, at radial distances beyond a few kpc. Hence, the 
most likely regime for the jets' decollimation due to the 
hot spots having slowed to subsonic speeds lies within 
roughly 10 kpc of the core. If the hot spot manages to 
retain supersonic speed out to such radial distances, then 
in most cases the jet is likely to propagate ahead down the 
rapidly declining ISM density, preserving its FR II char- 
acter. In some cases another possibility for jet flaring can 
arise farther out when the jet crosses the pressure-matched 
interface between the ISM and the intergala ctic m edium 
(GK W88; Wiita et al. [L990[ Wiita fc N orman |l992[ Hood a 
et al. |1994| ; Hooda & Wiita |l996j , |l998[ Zhang et al. |1999[ ). 
Relatively few sources will, however, make a transition to 
FR I morphology upon crossing this interface since a sig- 
nificant degree of beam collimation appears to occur by 
such radial distances within the extended halos of most 



galaxies (Blundell et al. 2001; Jeyakumar & Saikia 2000). 
For such very extended sources, our model is clearly over- 
simplified, as the evolution appears to break from self- 
similarity and the decrease in Pr with time (even for con- 
stant Lb) cannot be ignored (e.g., Blundell & Rawlings 

pooo| ). 

Other types of sources that may not be readily ac- 
commodated in our simple picture are, for instance, cases 
where a narrow jet is seen to terminate in a bright knot far 
away (say > 50 kpc) from the nucleus, and flares only be- 
yond that point (e.g., 3C 130, Hardcastle 1998). Of course, 
higher values of D* are expected for sources for which the 
beam's power and opening angle lie near the upper and 
lower ends of their respective ranges. Also, the possibility 
of a bright knot and jet expansion occurring when the jet 
interacts with irregularities in t he am bient medium must 
be considered (e.g., Wang et al. |2000 ). In this context it is 
especially interesting to recall that recent radio observa- 
tions of the nearest radio galaxy, Centaurus A, have pro- 
vided clear evidence that a radio jet can even emerge rec- 
ollim ated a fter flaring upon hitting a gas cloud (Morganti 
et al. 1999). These Australian Telescope Compact Array 



observations have revealed that the flaring of the north- 
ern jet of Cen A, a few kiloparsecs away from from the 
nucleus, gives rise to the inner radio 'lobe', out of which 
a narrow jet is seen to emerge and propagate further out 
for another few kiloparsecs before flaring and fading away 
in the fashion of typical FR I jets. Morganti et al. ( 1999[ ) 
have interpreted this striking 'leaky bulb' type morphol- 
ogy as being a manifestation of the de Laval nozzle model 
proposed by Blandford & Rees (1974). 

In our analytical formulation, we evaluate the critical 
beam power, L^, for which the hot-spot deceleration to 
subsonic velocities occurs at a fiducial distances, D*, of 3 
and 10 kpc from the core, for a range of absolute magni- 
tudes of the host galaxy between Mb = —19 and —23.5. 
The 10 kpc value is a typical distance at which jet s flar e 
in a sample of radio galaxies (O'Donoghue et al. 1993 ). 
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A recent study of a sample of 38 FR I sources by Laing 
et al. ( 1999Q gives a mean projected value of ~ 3.5 kpc 
for the radial distance of the point where the kiloparsec 
scale jet first becomes visible, after an initial emission gap 
(using our value of Hq). Our main results are not very 
sensitive to the choice of D* , at least until the very most 
luminous galaxies arc considered, and the predicted slope 
of the separation between FR I and FR II sources that we 
derive below is not strongly dependent upon it. 

We now write the relevant expressions as follows: 



X = C 2 n, 



2"0 



-1/2^1/2 



(6) 



where C 2 = (A/ncWrripiJ,) 1 ' 2 = 2.14 x 1CT 3 for fj, = 0.620. 
We also require, C\ = 47r(0.909)A(T)/x, where the cool- 
ing function, A(T) is a rather weak function of the temper- 
ature; assuming T = 10 7 K, A(T) ~ 1.17 x 10" 23 (Bicknell 



order approximation that the monochromatic radio lumi- 
nosity is a constant fraction of the beam power (at least 
near the transition luminosity). Then the local slope re- 
mains the same, with P£ 1 4GHz oc L^ p 5 t 6 , and the key fea- 
tures of the Owen-Ledlow diagram are satisfactorily ex- 
plained, with a very good fit to the slope retained. The 
intercept also fitted well, provided the product en ~ 0.27 
for D* = 10 kpc. As shown in Fig. 1, this is in excellent 
agreement with the bifurcation between FR I and FR II ra- 
dio sources, and nominally provides even a better fit than 
does Bicknell's (1995) model, which yields a power-law 
index above 2.0. We use a mean color index for elliptical 
galaxies (B — R) = 1.15 taken from a CCD photometric 



1995 ), and the integral over the gas density distribution 
gives Ix — 5.0 (for 6 = 0.75), a value we adopt, following 
Bicknell ( |1995| ). With these choices, d = 6.68 x 10" 22 . 
Equations (1-3) and (6) can now be combined to yield: 



log L x = log C\ + 21og n + 31og a, 



(7) 



logn x = logn -log«; = 4.188 + 0.218M B , (8) 

log X = -0.109M B +logC 2 +0.251ogCi+3.20+0.51ogL & .(9) 

We next impose the transonic condition on the LHS of Eq. 
(4), via 



v(D) =c s = [ 7 2.2/<5] 1 / 2 a. 



(10) 



Taking the ratio of specific heats for the combined plasma, 
7 = 5/3, and using 5 = 0.75, this becomes 

log v{D) = 0.345 + log a = -4.720 - 0.0959A/ B . (11) 

Equating this to the RHS of (4), treated as a function 
of X, gives, 



dex(-4.720 - 0.0959Af B ) 



X 



r, (12) 



D/[l + (D/a) 2 ] s / 2 +X' 
as the relation that must be solved for L£ in terms of M B . 

3. Results and Conclusions 

The values of L£, a, and no can be found for specified 
values of M B and D* by solving Eqs. (3), (8), (9) and 
(12). The results are given in Table 1, along with, (3, the 
local slope of the logL^-logL opt relation. For the choice of 
D* = 10 kpc, these slopes, (3, vary little for M B > -23.5, 



and correspond closely to a relation where Li oc L 



opt 



for 



D* = 3 kpc, (3 is only slightly less steep. For the most 
luminous galaxies, the value of from the model flattens 
somewhat. 

In the absence of a definitive approach to compute the 
efficiency, e, of conversion of the beam power into syn- 
chrotron radiation, we make the commonly adopted first 



study of cluster galaxies (J0rgensen et al. 1992) to convert 
from the absolute blue magnitudes used in Eqs. (1-3) to 
the red magnitudes plotted by Ledlow et al. (2001); the 



dependence of color index on optical magnitude is small 
over this range and we have not attempted to correct for 
it (J0rgensen et al. 1992; Kodama et al. 1998] ; Peletier et 
al. [19901 ) . 

We note that earlier estimates based on the application 
of beam models for FR II sources are in reasonable agree- 
ment with such values of e « 0.1 (Gopal-Krishna & Sari- 



palli 1984; Dreher 1984; Saripalli & Gopal-Krishna 1985), 



particularly if a "ready to radiate" electron-positron jet 
composition is assumed; and evidence f or such a compo- 
sition is g rowin g (e.g. 
Alexander 1997; 



Reyn olds e t al 
Kaiser et al 



1996a; Kaiser & 



1997; Hirotani et al. 1999) 



Given that our analysis involves a fairly large set of empir- 
ical relations for elliptical galaxies, the agreement of our 
prediction with the observational data is quite encourag- 
ing (Fig. 1). 

It is perhaps worth remarking on the consistency of our 
model with the recent Very Long Baseline Interferometry 
measurements for the expansion of very small and young 
Compact Symmetric Obj ect (CSO) radio sources (e.g., 
Owsianik & Conway gggg (OC)). For the CSO 0710+439, 
the best values obtained by OC are V — 0.13/i _1 c at a 
distance D ~ 25/i _1 pc (-C a), for uq ~ 2 cm" 3 , and 
L;, ~ 5 x 10 44 erg s _1 . Inserting the last three values into 
Eq. (4), and choosing our nominal 8 — 0.1 rad, we obtain 
V = 0.18h~ 1 c, in quite reasonable agreement; a somewhat 
higher value for 6 might be inferred from the maps in OC, 
thereby reducing our computed velocity and improving 
the agreement. Less well determined values of expansion 
speeds obtained by other groups for 3 other CSO radio 
sources quoted by OC range from 0.07 — 0.13/i _1 c and 
would easily be accommodated by more typical values of 
L b . The derived range, 0.08 < e < 0.31, for 0710+439 
(OC) is also in good agreement with our model. 

The forgoing analysis is admittedly oversimplified; for 
instance, we have considered a fixed opening angle for 
beams of all powers, and have also assumed specific values 
of 5 and k, all of which are certain to vary somewhat from 
source to source. On the other hand, the division between 
FR I and FR II sources is not perfect by any means, and 
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such a spread in properties would allow for the small num- 
ber of sources found on the 'wrong' side of the dividing 
line in Fig. 1. Also, we have ignored any variation in the 
jet disruption length with the galaxy's optical luminosity 
by typically choosing D* = 10 kpc; however, such a varia- 
tion is quite plausible, and therefore we have also provided 
results for D* = 3 kpc in the last two columns in Table 1. 

To obtain essentially as good a fit with D* = 3 kpc, 
one would have to raise e to 0.16 if k = 3 instead of 
e = 0.09 for D* — 10 kpc; of course both values could 
drop if k > 3 is considered. As a first step towards a more 
realistic model, one could consider a distribution of 2-10 
kpc for D* over the range in the hosts' optical magnitudes, 
taking smaller D* for lower L opt . Then the predicted slope 
of the Pft — L opt relation clearly would be slightly steeper 
(cf. columns 4 and 6 of Table 1), further improvi ng the 
accord with the data presented by Led low fc Owen ( 1996 ) 
and recently updated by Ledlow et al. ( 2001 ). We also note 
that the exact results quoted assume a specific spectral in- 
dex (a — 1.0) and upper and lower cut-offs to the radio 
band (10 MHz and 100 GHz) respectively, in relating the 
monochromatic radio power to the beam power. Assuming 
a flatter overall spectrum, say a = 0.7, would not affect 
the slope of the relation but would demand that the prod- 
uct en rise to ~ 0.52 (for D* = 10 kpc) to normalize it. 
Despite these uncertainties, it is probably fair to note that 
overall our variant model employs several fewer parame- 
ters than does Bicknell's (1995), although in connecting 



optical and X-ray properties we have closely followed his 
approach. 

The available very detailed maps of a few kiloparsec- 
scale FR I jets provide evidence for the existence of a 
slower moving sheath of synchrotr on plasma surround- 



ing a rel ativistic spine (e.g., Laing 1993 1996 ; Laing et 
al. 1999). In Bicknell's picture this sheath arises from 



turbulent mixing of the ambient plasma after the jet's 
Mach number becomes low enough to allow an exponen- 
tial growth of the planar Kelvin-Helmholtz instability. In 
our picture, the entrainment process is accelerated once 
the hotspot's motion has become subsonic, as this results 
in the cessation of the 'backflow' of the jet plasma; this 
leads to the diminution of the protective cocoon around 
the jet core which had hitherto separated the latter from 
the ISM material. We also recall that the impingement of 
the 'backflow' onto the jet engenders some of the recon- 
fmcment shocks (e.g. Norman et al. 1982; Hooda & Wiita 



1996, 1998) so that the depletion of the shrouding cocoon 



plasma reduces their effects and therefore contributes to 
the jet's decollimation. Thus, in our picture, the dramatic 
slowing down of the beam flow can be thought of as an 
eventual outcome of the decollimation of the beam's head 
(leading to FR I structure) , instead of being its principal 
cause. 
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Fig. 1. Radio flux plotted against host galaxy magnitude, 
with the numerals indicating FR class; data taken from 
Ledlow ct al. (2000). The superimposed curve is the di- 
vision between FR I and FR II sources arising from our 
model, taking D* = lOkpc and a constant fraction of to- 
tal beam power converted into synchrotron emission of 
e = 0.09. 



Table 1. Galaxy and jet parameters 
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